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Abstract 

In this note we formulate and investigate theoretical uncertainties for high Q 2 deep inelas- 
tic heavy quark (charm, etc.) production rates which arise within collinear resummation 
techniques from variations of the a priori unknown charm input scale Qo of 0(a s ) variable 
flavour number schemes. We show that Qo variations constitute a source of considerable 
theoretical uncertainty of present 0(a s ) calculations within such schemes and we suggest 
to consider a scale optimization from higher order corrections. We also outline how the 
stability of the fixed order and collinearly resummed perturbation series for heavy quark 
production can be comparatively investigated by variation of Qq. 



The present discussion on the appropriate scheme for the perturbative treatment of the 
deep inelastic production of heavy quarks of mass m 3> Aq C d can be partly traced back 
to the question what is the effective expansion parameter for high Q 2 predictions. While 
fixed order perturbation theory (FOPT) proceeds strictly stepwise in powers of (a s /27r) 
at all scales, variable flavor number schemes (VFNSs) are based upon the expectation 
that terms ~ (a s /2iT x InQ 2 /m 2 ) n from collinear regions in the phase space have to be 
resummed to all perturbative orders n for high Q 2 when (a s /2ir x \nQ 2 /m 2 ) — > 0(1). 
Such terms are undebatedly present in the high Q 2 limit of the perturbative partonic 
coefficient functions but their impact is less clear 0] on observable hadronic quantities 
like the charm component of the deep inelastic structure function F 2 C where the partonic 
coefficient functions have to be convoluted with modern, i.e. steep, parton distribution 
functions q (x,fi 2 F ) and dominantly g(x,fi 2 F ), fij ~ m. The question which ordering of 
the perturbation series optimizes its convergence can therefore not be answered a priori 
but only from explicit quantitative, i.e. numerical investigations 0,0]; prominent tools 
for testing perturbative stability being K factor considerations @ or scale variations 
||. At present both criteria indicate a well behaved fixed order perturbation series for 
relevant subasymptotic but large scales Q ^> m As regards scale uncertainties, 

mainly variations of the mass factorization scale /i^ have been considered so far despite 
the fact that collinear resummation techniques introduce an additional arbitrary scale in 
the process set by the input scale Q for the heavy quark: 



Recently proposed variable flavour number schemes p|-|TT| for global PDF analyses are 
constructed upon the boundary condition 



q { : f+1 \*,Qi: 







for a heavy sea quark density to enter the massless partonic renormalization group (RG)- 
evolution equations which resum collinear splitting subdiagrams to all orders at the price 
of neglecting mass dependent terms. In Eq. ([I]) m is the heavy quark mass and the heavy 
quark input scale Qo is in more technical terms the transition (or switching) scale from 
a factorization scheme with rif to the one with rtf + 1 partonic quark degrees of freedom 
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TT| . Since the scale Qq is of no physical meaning, a RG-like equation 

holds ideally for any heavy quark observable O. At limited perturbative order, Eq. (|2]) 
will obviously be violated to some extent which we will investigate below for the charm 
contribution to the NC structure function O = F£. 

At the heart of the variable flavour number schemes of PHTOH is some interpolation 
prescription between fixed order perturbation theory, assumed to be valid around Q 2 = 
0(m 2 ), and the Q 2 3> m 2 massless parton (MP) asymptotics derived from the boundary 
condition in Eq. ([I]). To avoid within our rather general considerations a discussion of 
the peculiarities of the distinct heavy quark schemes we denote such interpolations very 
schematically as 

{1, m 2 /Q 2 -> 1 
(3) 
0, m 2 /Q 2 -> 

where the simple weight w is meant to represent all the details of some elaborate scheme 
prescription P-ffD,[I^|. The deviation of VFNS from FOPT is thus normalized to MP 
and the predictive power of VFNS in Eq. @ depends on the stability of the asymptotic 
MP prediction which is obtained from the boundary ([!]) at Qq = m via massless RG 
evolutions. Equation (|I]) emerges from the matching conditions of a factorization scheme 
with n j active flavours to a scheme with nj + 1 active flavours at some a priori arbitrary 
transition (or switching) scale Qq. The general transformation equations for quark (q) 
and gluon (g) parton densities as well as for a s read up to NLO jO], |l^] 

1 { : f+1) M 2 a) = ^ m§ £ f P$>(0 (f, Q§) + 0(a 2 ) 
9^ +1) (x,Q 2 ) = g^(x,Q 2 ) (l + -Ml Ing) + 0{a 2 s ) 

a^ +1) (Ql) = a^\Q 2 ) /(l + ^lnf) + O(o2) 

q^ +1) (x,Q 2 ) = q^\x,Q 2 ) + 0(a 2 ) 



(4) 



and obviously reduce to Eq. ([]]) for Q 



m: 



Q 

1 n 



(nf+1 \x,m 2 ) = , g (n f +1 \x,m 2 ) = g^\x,m 2 ) , af s+l \m 2 ) = a { : f \m 2 ) . (5) 



In Eqs. (^) and (pD q introduces a partonic heavy quark density into the massless evolu- 
tion equations and the unspecified a s (Ql) may be either a^^iQo) or «i n/+1 ^(Qo) because 
the difference is of the orders neglected in (f|). The argument of a continuous g(x, Ql) 
and a s {Qo) has been advanced 0,|lTJ for adopting Q = m in all NLO parton distribu- 



tion sets constructed so far along a VFNS philosophy |T3|-|T3|. On the other hand, when 
Qo 7^ m the discontinuity of a s is in practice as small as maximally 4% up to Q$ as high 
as 1000 GeV 2 . Anyway, the recently completed NNLO transformation equations 
reveal that the possibility of a continuous evolution across Qo breaks down beyond NLO 
by nonlogarithmic higher order corrections to (|4]), (|5|). The restriction to Qo = m should 
hence be abandoned and effects of varying Q Q should be taken into account on the same 
level as the variations of the mass factorization scale \j? F usually considered. Indeed, along 



ln^ = In ^ + In 4 + In % (6) 

m 2 m 2 Qq \ip 



the two scales Qo and \ip define quite symmetrically which portion of the quasi-collinear 
\\\{Q 2 /m 2 ) is actually resummed [\a(fiF / Qo)] and what amount is kept at fixed order, 
either in the boundary condition for q n [ln(Ql/m 2 )] or in the hard scattering coefficient 
function Cf in Eq. (0) below [ln(Q 2 /^)]. We will investigate the residual Qo dependence 
for the charm production contribution to the deep inelastic structure function F<i using 
m = m c (= 1.5 GeV). To avoid complications from an interplay of several scales we will 
decouple the bottom and top quark from the process (m^ — > oo) and we will fix the 
factorization scale at Hf — Q- 



In the asymptotic limit m 2 c /Q 2 — > the schemes [0,|9|,[T(J reduce - as in Eq. (|3|) - 
to the so-called Zero Mass Variable Flavour Number Scheme, equivalent to the 'massless 
parton' scenario of Ref. where any mass dependence is dropped except for the boundary 
conditions in (f|). We will consider such a scenario in the following and ignore terms of 
0{m 2 /Q 2 ) because these are not handled uniformly in the individual realizations J7|,P|.pH 
of a VFNS.Q For definiteness we consider an F^ MP in 7*P scattering which is given by 

^ MP (x,Q 2 ) = e 2 c \(c + c)(x,Q 2 ) + ^4r^{9®C! + (c + c^ , (7) 

X Z7T 



indeed the charm scheme of JTo[ is explicitly constructed upon Qo = m and a generalization of this 
particular scheme for Qo ^ rn seems nontrivial. 
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Figure 1: Dependence of the 'massless parton'-type charm structure function _F 2 C,MP in 
Eq. (|7|) on the switching scale Qq in the boundary conditions in Eq. (|4]). The solid lines 
represent the central value of Qq = m c and F 2 C,MP decreases monotonically with increasing 



Qq. The underlying parton distributions below Q are those of Ref. |13| . 



where the Cf' 9 are the massless MS coefficient functions (T^,[T8[. It has already been 



pointed out in || that the 'massless parton' F2 MF in Eq. (^) can be rather arbitrarily 



suppressed if some larger effective charm mass is introduced |19| into the boundary con- 
dition (|^). Our investigation here will clarify the situation if - for a fixed value of the 
physical charm mass m c - the unphysical switching scale Qq is varied consistently accord- 
ing to the NLO boundary equations (f|). Fig. [I] shows the effect if the transition scale 
is allowed to vary over the range m 2 c /2 < Q^ < 2 m? c where the central value Qq = m c 



is represented by the solid lines and _F 2 C ' monotonically decreases with increasing Qo-Q 
The evolution leading to the results in Fig. [I] is based on the rif = 3 valence-like NLO 



input of Ref. ]T3[ using NLO (2-loop) splitting functions. Above Qq the evolution deviates 
from |l3j because we consider general Qq ^ m c here and we ignore - as mentioned above - 



any bottom quark effects (m^ — > oo). The amount of change of F^ MP under variation of 
Qq hints at a reasonable perturbative stability. Nevertheless, the error represented by the 



2 Allowing for Qo < m c in Eq. (||) leads obviously to c(x, Qq) < which appears somewhat counter- 
intuitive in probabilistic parton model language. Note, however, that a negative charm input arises even 
for Qo = m from higher order corrections to Eq. (|j) |2Cf| . Anyway, the measurable cross section F% is 
certainly positive above the physical threshold Q 2 (l/x — 1) > 4m 2 . 
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Figure 2: Dependence of the 'massless parton'-type charm structure function _F 2 C,MP in 
Eq. (0) on the switching scale Q in the boundary conditions in Eq. (f|). In this Figure Q 
is allowed to vary (maximally) between the input scale of pi as a lower and the physical 
scale Q 2 as an upper limit. 



shaded bands in Fig. [I] is of the typical order of discrepancies between VFNS and FOPT 



calculations |10|,|T2[ which questions the gain in predictivity if FOPT is abandoned for 
VFNS. Such uncertainties are critical for precise charm predictions to compare with fu- 
ture experimental accuracy, especially at experimentally most relevant intermediate scales 
and regarding the fact that Qq was only allowed to fluctuate by a factor of 2. This latter 
limitation rests on the assumption that Qo has to be very close to m c for the all-order 
logarithms (a s /2ir x In Q 2 / 'm 2 c ) n to be correctly resummed. One can as well adopt a very 
different point of view towards the choice of Qq. One can easily see that inserting Eq. (^) 
into Eq. (0) gives 



lim F- MP (x,Q 2 ) 

Qo—>Q 



?GF 



1 + T^ 



[X 



,Q 2 ) + 0(a 2 s ) + 0[^ 



nr 



(8) 



which is dominated by the 0[a s \n(Q 2 /m 2 )] LO gluon fusion term F^ F of the fixed order 
perturbation series. We may hence consider Qo — > Q in a sense as a continuous path from 
variable flavour number to fixed order calculations. We should then consider values of 
Qo as high as we trust fixed order perturbation theory. In Fig. |2| we cover the maximally 
conceivable range for Qq] i.e. the leftmost end of all curves is set by the low input scale of 
the parton distributions in [O], Ql = 0.3 GeV 2 , while the rightmost ends are the 'fixed 



order limit' in Eq. (|8|). The observed monotonic scale dependence has to be expected from 
the positivity of the collinear resummation at small x which is continuously suppressed the 
more Qq is increased. Worrisome is, however, the steep slope dF2 MP /d\n.{Ql/ra 2 c ) around 
Qq ~ m c for high W 2 = Q 2 (l/x — 1), where the charm contribution is most important. 
This observation restricts the predictive power of NLO collinear resummation techniques 
which have thus far been constructed to match the asymptotic (Q 2 — > oo) F^ MP derived 
from c(x, Qq = m 2 ) = 0. The uncertainty from the residual Qq dependence, inherent to 
any VFNS ||7|, |9|, [10|] worked out to NLO, seems to dominate over scheme uncertainties of 
0(m 2 /Q 2 ) between the individual schemes. The arbitrariness of Qq therefore constitutes 
a main limiting factor on the perturbative accuracy of VFNS heavy quark predictions at 
high Q 2 . On the other hand we observe a flattening slope <9.F 2 C ' MP /d\\i{Ql/m 2 c ) towards 
the 'fixed order limit' at high Qq < Q where perturbative NLO <-> LO stability had been 
found in by K factor considerations for the full fixed order predictions, i.e. including 
logs and finite terms. We should reemphasize that these conclusions are based on the NLO 
matching conditions in Eq. (f|) - which is the present state of the art for PDF sets including 
partonic heavy quarks |T3|-[I5| - and do not take into account the higher corrections of 
12|,[L6|]. The terms beyond Eq. (ffl) represent NNLO contributions to the asymptotic 



VFNS prediction F^ MP ■ Very recently the results of [12] have been implemented in a 
0(a 2 ) implementation of a VFNS0 where the contribution from the unknown NNLO (3- 
loop) splitting functions had to be neglected. Choosing Qq = m the results of [^(J seem 
to indicate that the impact of terms from the resummation beyond fixed NLO [0(a^)]F] 
perturbation theory is rather moderate. As a further step in the line of the present 



investigation it would clearly be interesting to generalize p0[ to Qq ^ m. If a Qq > m 
would be prefered by such an analysis the difference between VFNS and FOPT would be 
reduced even more. Such a result would again re-confirm the perturbative reliability of 
FOPT found in as much as it would help reduce unphysical scheme dependences of 
QCD predictions on charm production and thus make a comparison to experiment even 



3 A c(x,Q 2 ) derived from the NNLO boundary conditions in fT^I would, however, be problematic to 
apply to hadroproduction calculations, since the higher terms in |[L2| are not yet contained in, e.g., the 
fixed NLO [0(o!g)] hadroproduction process pp — > c{pt)X |f2l| . 

The confusion of counting perturbative orders differently in resummed (MP) and fixed order (FOPT) 
calculations is treated in more detail in H, [10| . 
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more compelling. 



To summarize, we have considered variations of the a priori arbitrary charm input scale 
Qo, which separates a 3 from a 4 flavour scheme in variable flavour number approaches, 
around its usually adopted but by no means theoretically required value of Qq = m c . From 
the NLO boundary conditions we found a monotonic Qo dependence and a worrisome steep 
slope of F2 MP in (0) with respect to lnQp just around Qo ~ m c- This behaviour restricts 
the accuracy of a collinearly resummed NLO approach towards calculating high Q 2 charm 
production from matching 0(a s ) boson gluon fusion at Qq = m c in the MP component of 
Eq. (||). This uncertainty in MP from the unknown Q feeds back onto the entire VFNS 
via Eq. (^) which is normalized to MP asymptotically (Q 2 ^> m 2 ). Our results imply to 
consider variations of Qo both as a limiting factor on the present perturbative accuracy if 
estimating the theoretical uncertainty of VFNS heavy quark predictions as well as in order 
to optimize the starting scale for the charm evolution within higher order realizations of 
VFNSs P"2| , P0| , . This latter higher order analysis assumes, however, that the unknown 
NNLO (3-loop) splitting functions can be neglected |20 . 
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